2. Kinetic Theory 



The purpose of this section is to lay down the foundations of kinetic theory, starting 
from the Hamiltonian description of 10^^ particles, and ending with the Navier-Stokes 
equation of fluid dynamics. Our main tool in this task will be the Boltzmann equation. 
This will allow us to provide derivations of the transport properties that we sketched in 
the previous section, but without the more egregious inconsistencies that crept into our 
previous derivaion. But, perhaps more importantly, the Boltzmann equation will also 
shed light on the deep issue of how irreversibility arises from time-reversible classical 
mechanics. 

2.1 From Liouville to BBGKY 

Our starting point is simply the Hamiltonian dynamics for N identical point particles. 
Of course, as usual in statistical mechanics, here is N ridiculously large: N ~ 0(10^^) 
or something similar. We will take the Hamiltonian to be of the form 

N N 

^ - ^ E^t + E ^(^^) + E ^(^"^ - (2-1) 

i=l i=l i<j 
— * 

The Hamiltonian contains an external force F — —W that acts equally on all parti- 
cles. There are also two-body interactions between particles, captured by the potential 
energy [/(r^ — fj). At some point in our analysis (around Section 2.2.3) we will need 
to assume that this potential is short-ranged, meaning that C/(r) ~ for r ^ d where, 
as in the last Section, d is the atomic distance scale. 

Hamilton's equations are 

dpi_ ^ _dH_ dn^dH_ 
dt dfi ' dt dpi 

Our interest in this section will be in the evolution of a probability distribution, 
f{fi,pi;t) over the 2N dimensional phase space. This function tells us the proba- 
bility that the system will be found in the vicinity of the point {fi,pi). As with all 
probabilities, the function is normalized as 

N 

/ dVf{fi,pi; t) = 1 with dV = l[ d\d^pi 

i=l 

Furthermore, because probability is locally conserved, it must obey a continuity equa- 
tion: any change of probability in one part of phase space must be compensated by 
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a flow into neighbouring regions. But now we're thinking in terms of phase space, 
the "V" term in the continuity equation includes both d/dvi and d/dpi and, corre- 



This flnal equation is the Liouville's equation. It is the statement that probabihty 
doesn't change as you follow it along any trajectory in phase space, as is seen by 
writing the Liouville equation as a total derivative. 



To get a feel for how probability distributions evolve, one often evokes the closely 
related Liouville's theorem^. This is the statement that if you follow some region of 
phase space under Hamiltonian evolution, then its shape can change but its volume 
remains the same. This means that the probability distribution on phase space acts 
like an incompressible fluid. Suppose, for example, that it's a constant, /, over some 
region of phase space and zero everywhere else. Then the distribution can't spread out 
over a larger volume, lowering its value. Instead, it must always be / over some region 
of phase space. The shape and position of this region can change, but not its volume. 

The Liouville equation is often written using the Poisson bracket, 



spondingly, the velocity vector in phase space is (fijPj). The continuity equation of the 
probability distribution is then 




But, using Hamilton's equations (2.2), this becomes 





OA dB dA dB 

dfi dpi dpi dfi 



With this notation, Liouville's equation becomes simply 



df 



{HJ} 



dt 



fuller discussion of Hamiltonian mechanics and Liouville's theorem can be found in Section 4 of 
the classical dynamics notes: http://www.damtp.cam.ac.uk/user/tong/dynamics.html. 



It's worth making a few simple comments about these probabihty distributions. Firstly, 
an equilibrium distribution is one which has no explicit time dependence: 



which holds if {H, /} = 0. One way to satisfy this is if / is a function of H and the most 
famous example is the Boltzmann distribution, / ~ e~^^ . However, notice that there 
is nothing (so-far!) within the Hamiltonian framework that requires the equilibrium 
distribution to be Boltzmann: any function that Poisson commutes with H will do the 
job. We'll come back to this point in Section 2.2.2. 

Suppose that we have some function, A{fi,pi), on phase space. The expectation 
value of this function is given by 



This expectation value changes with time only if there is explicit time dependence in 
the distribution. (For example, this means that in equilibrium (A) is constant). We 
have 



where we have integrated by parts to get to the last line, throwing away boundary 
terms which is justified in this context because / is normalized which ensures that we 
must have / — > in asymptotic parts of phase space. Finally, we learn that 



This should be ringing some bells. The Poisson bracket notation makes these expres- 
sions for classical expectation values look very similar to quantum expectation values. 

2.1.1 The BBGKY Hierarchy 

Although we're admitting some ignorance in our description of the system by consider- 
ing a probability distribution over iV-particle phase space, this hasn't really made our 
life any easier: we still have a function of ~ 10^^ variables. To proceed, the plan is 




(2.3) 




(2.4) 




(2.5) 
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to limit our ambition. We'll focus not on the probability distribution for all N parti- 
cles but instead on the one-particle distribution function. This captures the expected 
number of parting lying at some point (r , p). It is defined by 



N 



/.(r>;0^7v/n<^VA/(r-f.....r-„.p>-....p.;0 

Although we seem to have singled out the first particle for special treatment in the 
above expression, this isn't really the case since all N of our particles are identical. 
This is also refiected in the factor N which sits out front which ensures that /i is 
normalized as 



d-'rdy h{rj-t) = N (2.6) 

For many purposes, the function /i is all we really need to know about a system. In 
particular, it captures many of the properties that we met in the previous chapter. For 
example, the average density of particles in real space is simply 



n{f-t) = J d'ph{fj-t) (2.7) 
The average velocity of particles is 

u{r-t)^ I d^p^Mr,p;t) (2.8) 

and the energy flux is 

E{r; t)^ [ d'p ^E{p)Mf,p; t) (2.9) 
J m 

where we usually take E{p) = p^ /2m. All of these quantities (or at least close relations) 
will be discussed in some detail in Section 2.4. 

Ideally we'd like to derive an equation governing /i. To see how it changes with time, 
we can simply calculate: 

^ = AT 1 n d"^^d^P^ d'nd'P^ {H, /} 

Using the Hamiltonian given in (2.1), this becomes 



dt 



„ N 

- = AT / l[d\d^p, 

1=2 



^ m dfi ^ dfi dp - ^ ^ 



^ j=i k<i ^ 
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Now, whenever j = 2, . . . N, we can always integrate by parts to move the derivatives 
away from / and onto the other terms. And, in each case, the result is simply zero 
because when the derivative is with respect to fj, the other terms depend only on pj 
and vice-versa. We're left only with the terms that involve derivatives with respect 
to fi and pi because we can't integrate these by parts. Let's revert to our previous 
notation and call f i = r and pi = p. We have 



^-Njfld^.d^p 



p df ^ dV{f) df ^ dU{f- fk) df 



m df df dp df dp 

= {HuM + Nff[ ^r.Sp. ± ■ I (2.10) 

i=2 k=2 ^ 

where we have defined the one-particle Hamiltonian 

Notice that Hi includes the external force V acting on the particle, but it knows nothing 
about the interaction with the other particles. All of that information is included in 
the last term with U{f—fk). We see that the evolution of the one-particle distribution 
function is described by a Liouville-like equation, together with an extra term. We 
write 

df 
dt 



The first term is sometimes referred to as the streaming term. It tells you how the 
particles move in the absence of collisions. The second term, known as the collision 
integral, is given by the second term in (2.10). In fact, because all particles are the 
same, each of the {N — 1) terms in '^^=2 (2-10) are identical and we can write 

(^)co,r "^^"^ " / ^^^^w^ ■ Yp J n^'-A/(r,f., . . . . . . ;t) 

But now we've got something of a problem. The collision integral can't be expressed 
in terms of the one-particle distribution function. And that's not really surprising. As 
the name suggests, the collision integral captures the interactions - or collisions - of 
one particle with another. Yet /i contains no information about where any of the other 
particles are in relation to the first. However some of that information is contained in 
the two-particle distribution function, 

N 

f2{fi,r2,Pi,P2;t) = N{N - 1) / Yld^rid^pi f{fi,f2, . . . ,pi,P2, . . .■,t) 

1=3 
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With this definition, the coUision integral is written simply as 

The coUision term doesn't change the distribution of particles in space. This is captured 
by the particle density (2.7) which we get by simply integrating n — J d^pfi- But, after 
integrating over J d^p, we can perform an integrating by parts in the collision integral 
to see that it vanishes. In contrast, if we're interested in the distribution of velocities - 
such as the current (2.8) or energy fiux (2.9) - then the collision integral is important. 

The upshot of all of this is that if we want to know how the one-particle distribution 
function evolves, we also need to know something about the two-particle distribution 
function. But we can always figure out how /2 evolves by repeating the same calculation 
that we did above for /i. It's not hard to show that /2 evolves by a Liouville-like 
equation, but with a corrected term that depends on the three-particle distribution 
function f^. And evolves in a Liouville manner, but with a correction term that 
depends on f^, and so on. In general, the n-particle distribution function 

ATI f ^ 

obeys the equation 

-m = + 1^ y ^ '"^^^ ^"^^ dr. W ^'-'^^ 

where the effective n-body Hamiltonian includes the external force and any interactions 
between the n particles but neglects interactions with any particles outside of this set. 



The equations (2.14) are known as the BBGKY hierarchy. (The initials stand for 
Bogohubov, Born, Green, Kirkwood and Yuan). They are teUing us that any group 
of n particles evolves in a Hamiltonian fashion, corrected by interactions with one of 
the particles outside that group. At first glance, it means that there's no free lunch; 
if we want to understand everything in detail, then we're going to have to calculate 
everything. We started with the Liouville equation governing a complicated function 
/ of iV ~ 0(10^^) variables and it looks like all we've done is replace it with 0(10^^) 
coupled equations. 
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However, there is an advantage is working with the hierarchy of equations (2.14) 
because they isolate the interesting, simple variables, namely /i and other lower /„. 
This means that the equations are in a form that is ripe to start implementing various 
approximations. Given a particular problem, we can decide which terms are important 
and, ideally, which terms are so small that they can be ignored, truncating the hierarchy 
to something manageable. Exactly how you do this depends on the problem at hand. 
Here we explain the simplest, and most useful, of these truncations: the Boltzmann 
equation. 

2.2 The Boltzmann Equation 

In this section, we explain how to write down a closed equation for /i alone. This 
will be the famous Boltzmann equation. The main idea that we will use is that there 
are two time scales in the problem. One is the time between collisions, r, known as 
the scattering time or relaxation time. The second is the collision time, Tcou, which 
is roughly the time it takes for the process of collision between particles to occur. In 
situations where 



we should expect that, for much of the time, /i simply follows its Hamiltonian evolution 
with occasional perturbations by the collisions. This, for example, is what happens for 
the dilute gas. And this is the regime we will work in from now on. 

At this stage, there is a right way and a less-right way to proceed. The right way is 
to derive the Boltzmann equation starting from the BBGKY hierarchy. And we will do 
this in Section 2.2.3. However, as we shall see, it's a little fiddly. So instead we'll start 
by taking the less-right option which has the advantage of getting the same answer 
but in a much easier fashion. This option is to simply guess what form the Boltzmann 
equation has to take. 

2.2.1 Motivating the Boltzmann Equation 

We've already caught our first glimpse of the Boltzmann equation in (2.12), 



But, of course, we don't yet have an expression for the collision integral in terms of 
/i. It's clear from the definition (2.13) that the second term represents the change in 
momenta due to two-particle scattering. When r ^ Tcou, the collisions occur occasion- 
ally, but abruptly. The collision integral should reflect the rate at which these collisions 
occur. 



T > Tcoll 



(2.15) 




(2.16) 
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Suppose that our particle sits at (r, in phase space and colhdes with another 
particle at {f,p2). Note that we're assuming here that collisions are local in space so 
that the two particles sit at the same point. These particles can collide and emerge 
with momenta p[ and We'll define the rate for this process to occur to be 

Rate = uj{p,P2\p[,P2) /2(r, r,p,P2) d^P2d^p'id^P2 (2.17) 

(Here we've dropped the explicit t dependence of /2 only to keep the notation down). 

The scattering function ou contains the information about the dynamics of the process. 
It looks as if this is a new quantity which we've introduced into the game. But, using 
standard classical mechanics techniques, one can compute uj for a given inter-atomic 
potential U{r). (It is related to the differential cross-section; we will explain how to do 
this when we do things better in Section 2.2.3). For now, note that lj is proportional to 
the two-body distribution function /2 since this tells us the chance that two particles 
originally sit in {f,p) and (r,p2)- 

We'd like to focus on the distribution of particles with some specified momentum 
p. Two particles with momenta p and p2 can be transformed in two particles with 
momenta p[ and Since both momenta and energy are conserved in the collision, we 
have 

P + P2=pi+P2 (2.18) 
p'+pI=p'i+p'2 (2.19) 

There is actually an assumption that is hiding in these equations. In general, we're 
considering particles in an external potential V. This provides a force on the particles 
which, in principle, could mean that the momentum and kinetic energy of the particles 
is not the same before and after the collision. To eliminate this possibility, we will 
assume that the potential only varies appreciably over macroscopic distance scales, so 
that it can be neglected on the scale of atomic collisions. This, of course, is entirely 
reasonable for most external potentials such as gravity or electric fields. Then (2.18) 
and (2.19) continue to hold. 

While collisions can deflect particles out of a state with momentum p and into a 
different momentum, they can also deflect particles into a state with momentum p. 
This suggests that the coUision integral should contain two terms, 

^ {p[ , P2 \P, Pi ) h {r, f,p[,p2)-u; {p, p2 \p[ , P2) /2 {r, f, p, P2) 

The first term captures scattering into the state p, the second scattering out of the 
state p. 



^ ) = / d'p2d'pW2 



mil 



The scattering function obeys a few simple requirements. Firstly, it is only non- 
vanishing for scattering events that obey the conservation of momentum (2.18) and 
energy (2.19). Moreover, the discrete symmetries of spacetime also give us some im- 
portant information. Under time reversal, p —p and, of course, what was coming in 
is now going out. This means that any scattering which is invariant under time reversal 
(which is more or less anything of interest) must obey 

^^{P,P2\P[,P2) = I^{-P[, -P2I -P, -P2) 

Furthermore, under parity {f,p) — > (— r, —p). So any scattering process which is parity 
invariant further obeys 

^{P,P2\P[,P2) = ^(-P, -P2I -P[, -P2) 

The combination of these two means that the scattering rate is invariant under exchange 
of ingoing and outgoing momenta, 

^{P,P2\p'l,P2) = ^{P[,P2\P:P2) (2.20) 

(There is actually a further assumption of translational invariance here, since the scat- 
tering rate at position —r should be equivalent to the scattering rate at position +r). 

The symmetry property (2.20) allows us to simplify the collision integral to 



d p2d p[d P2 u}{p[,P2\p,P2) f2{r,f,p[,P2) - f2{r,r,p,P2) (2.21) 



To finish the derivation, wc need to face up to our main goal of expressing the collision 
integral in terms of /i rather than /2. We make the assumption that the velocities of 
two particles are uncorrelated, so that we can write 

f2{r,r,P,P2) = fi{r,p)fi{r,P2) (2.22) 

This assumption, which sometimes goes by the name of molecular chaos, seems innocu- 
ous enough. But actually it is far from innocent! To see why, let's look more closely 
at what we've actually assumed. Looking at (2.21), we can see that we have taken the 
rate of collisions to be proportional to /2(^, P2) where pi and p2 are the momenta 
of the particles before the collision. That means that if we substitute (2.22) into (2.21), 
we are really assuming that the velocities are uncorrelated before the collision. And 
that sounds quite reasonable: you could imagine that during the collision process, the 
velocities between two particles become correlated. But there is then a long time, r, 
before one of these particles undergoes another collision. Moreover, this next collision 
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is typically with a completely different particle and it seems entirely plausible that the 
velocity of this new particle has nothing to do with the velocity of the first. Nonethe- 
less, the fact that we've assumed that velocities are uncorrected before the collision 
rather than after has, rather slyly, introduced an arrow of time into the game. And 
this has dramatic implications which we will see in Section 2.3 where we derive the 
H-theorem. 

Finally, we may write down a closed expression for the evolution of the one-particle 
distribution function given by 

with the collision integral 

(^) /^'^^2ciVl^iV2^(p^,P2|P,P2)[/l(r,pO/l(^ (2-24) 

This is the Boltzmann equation. It's not an easy equation to solve! It's a differential 
equation on the left, an integral on the right, and non-linear. You may not be surprised 
to hear that exact solutions are not that easy to come by. We'll see what we can do. 

2.2.2 Equilibrium and Detailed Balance 

Let's start our exploration of the Boltzmann equation by revisiting the question of the 
equilibrium distribution obeying df^'^/dt = 0. We already know that {f,Hi} = if / 
is given by any function of the energy or, indeed any function that Poisson commutes 
with H. For clarity, let's restrict to the case with vanishing external force, so V{r) = 0. 
Then, if we look at the Liouville equation alone, any function of momentum is an 
equilibrium distribution. But what about the contribution from the collision integral? 

One obvious way to make the collision integral vanish is to find a distribution which 
obeys the detailed balance condition, 

/r(r, Pl)/r(r, p1) = /r(r, p1/r(r, P2) (2.25) 

In fact, it's more useful to write this as 

log(/r(f,p1)) + log(/r(f,p1)) = log(/r(r,p-)) + log(/r(r,P2)) (2.26) 

How can we ensure that this is true for all momenta? The momenta on the right are 
those before the collision; on the left they are those after the collision. From the form of 
(2.26), it's clear that the sum of logf^"^ must be the same before and after the coUision: 
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in other words, this sum must be conserved during the colhsion. But we know what 
things are conserved during colhsions: momentum and energy as shown in (2.18) and 
(2.19) respectively. This means that we should take 



log(/r(r,p-)) ^(3{^x-E{p)+u■ p) (2.27) 

where E{p) = /2m for non-relativistic particles and /x, j3 and u are all constants. 
We'll adjust the constant /i to ensure that the overall normalization of /i obeys (2.6). 
Then, writing p = mv, we have 

/r(r, P) = y (^) e-^-(— (2.28) 

which reproduces the Maxwell-Boltzmann distribution if we identify /3 with the inverse 
temperature. Here u allows for the possibihty of an overall drift velocity. We learn 
that the addition of the coUision term to the Liouville equation forces us to sit in the 
Boltzmann distribution at equilibrium. 

There is a comment to make here that will play an important role in Section 2.4. 
If we forget about the streaming term then there is a much larger class of 

solutions to the requirement of detailed balance (2.25). These solutions are again of 
the form (2.27), but now with the constants /i, (5 and u promoted to functions of space 
and time. In other words, we can have 

fr\r,m = n(f,t) (§^)'^'exp (-/9(r,t)y [{v - u{r,t)f) (2.29) 

Such a distribution is not quite an equilibrium distribution, for while the collision 
integral in (2.23) vanishes, the streaming term does not. Nonetheless, distributions 
of this kind will prove to be important in what follows. They arc said to be in local 
equilibrium, with the particle density, temperature and drift velocity varying over space. 

The Quantum Boltzmann Equation 

Our discussion above was entirely for classical particles and this will continue to be 
our focus for the remainder of this section. However, as a small aside let's look at how 
things change for quantum particles. We'll keep the assumption of molecular chaos, 
so /2 ~ /i/i as in (2.22). The main difference occurs in the scattering rate (2.17) for 
scattering Pi + P2 ^ Pi + P2 which now becomes 

Rate - uj{pj2\p[J2) h{Pi)h{m 1 ± fM)} {1 ± /(P2)} d'p2d'p[d% 
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The extra terms are in curly brackets. We pick the + sign for bosons and the — sign 
for fermions. The interpretation is particularly clear for fermions, where the number of 
particles in a given state can't exceed one. Now it's not enough to know the probability 
that the final state is free for the particle to scatter into: and that's what the {1 — /i} 
factors are teUing us. 

The remaining arguments go forward as before, resulting in the quantum Boltzmann 
equation 

^) = [d'p2d'p[d'p',u;{p[,p^\p,p2) \hip[)Mp^){l±Mp)}{l±Mp2} 



coll 



-Mp)MP2) {1 ± {1 ± /i(r,p1)} 



To make contact with what we know, we can look again at the requirement for equi- 
librium. The condition of detailed balance now becomes 



li/rvo; "Vii/rvay "vii/r®; "vii/n??.). 

Which is again solved by relating each log to a linear combination of the energy and 
momentum. We find 

1 



which reproduces the Bose-Einstein and Fermi-Dirac distributions. 
2.2.3 A Better Derivation 

In Section (2.2.1), we derived an expression for the collision integral (2.24) using in- 
tuition for the scattering processes at play. But, of course, we have a mathematical 
expression for the collision integral in (2.13) involving the two-particle distribution 
function /2. In this section we will sketch how one can derive (2.24) from (2.13). This 
will help clarify some of the approximations that we need to use. At the same time, 
we will also review some basic classical mechanics that connects the scattering rate uj 
to the inter-particle potential U (r) . 

We start by returning to the BBGKY hierarchy of equations. For simplicity, we'll 
turn off the external potential V{r) = 0. We don't lose very much in doing this because 
most of the interesting physics is concerned with the scattering of atoms off each other. 
The first two equations in the hierarchy are 

^ + — ■ ^ /i = d r2d p2 ^ jpr 2.30 

at m ori ) j ori opi 
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and 

d d 



d ^Pi d ^ p2 d 1 dU {fi — 



dt m dfi m df2 2 Of 



dpi dp 



f2 = (2.31) 



dU{fi - fs) (9 dU{r2 - r^) ^ \ j 
dfi dpi df2 dp 



In both of these equations, we've gathered the streaming terms on the left, leaving 
only the higher distribution function on the right. To keep things clean, we've sup- 
pressed the arguments of the distribution functions: they are /i = fi{fi,pi;t) and 
/2 = f2{ri,r2,Pi,P2]t) and you can guess the arguments for f^. 

Our goal is to better understand the collision integral on the right-hand-side of (2.30). 
It seems reasonable to assume that when particles are far-separated, their distribution 
functions are uncorrelated. Here, "far separated" means that the distance between 
them is much farther than the atomic distance scale d over which the potential U (r) 
extends. We expect 

f2{ri,r2,Pi,P2;t) fi{fi,Pi;t)fi{f2,P2;t) when |fi -fi| > d 

But, a glance at the right-hand-side of (2.30) tells us that this isn't the regime of 

interest. Instead, /2 is integrated dU{r)/dr which varies significantly only over a region 
r < d. This means that we need to understand /2 when two particles get close to each 
other. 

We'll start by getting a feel for the order of magnitude of various terms in the 
hierarchy of equations. Dimensionally, each term in brackets in (2.30) and (2.31) is an 
inverse time scale. The terms involving the inter-atomic potential U (r) are associated 
to the collision time Tcou- 

1 dU d 
Tcou df dp 

This is the time taken for a particle to cross the distance over which the potential U (r) 
varies which, for short range potentials, is comparable to the atomic distance scale, d, 
itself and 

d 

Vrel 

where v^ei is the average relative speed between atoms. Our first approximation will be 
that this is the shortest time scale in the problem. This means that the terms involving 
dU/dr are typically the largest terms in the equations above and determine how fast 
the distribution functions change. 
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With this in mind, we note that the equation for fi is special because it is the only 
one which does not include any collision terms on the left of the equation (i.e. in 
the Hamiltonian This means that the collision integral on the right-hand side 

of (2.30) will usually dominate the rate of change of /i. (Note, however, we'll meet 
some important exceptions to this statement in Section 2.4). In contrast, the equation 
that governs /2 has collision terms on the both the left and the right-hand sides. But, 
importantly, for dilute gases, the term on the right is much smaller than the term on 
the left. To see why this is, we need to compare the term to the /2 term. If we were 
to integrate /a over all space, we get 

(where we've replaced {N — 2) Ri N in the above expression). However, the right- 
hand side of (2.31) is not integrated over all of space. Instead, it picks up a non-zero 
contribution over an atomic scale ~ c?^. This means that the collision term on the 
right- hand- side of (2.31) is suppressed compared to the one on the left by a factor of 
Nd^/V where V is the volume of space. For gases that we live and breath every day, 
Nd^/V ~ 10~^ — 10~^. We make use of this small number to truncate the hierarchy of 
equations and replace (2.31) with 

f d pi d P2 d 1 dU{fi — f2) 
\dt m dfi m df2 2 dfi 

This tells us that /2 typically varies on a time scale of Tcqu and a length scale of d. 
Meanwhile, the variations of fi is governed by the right-hand-side of (2.30) which, by 
the same arguments that we just made, are smaller than the variations of /2 by a factor 
of Nd^/V. In other words, fi varies on the larger time scale r. 

In fact, we can be a little more careful when we say that /2 varies on a time scale 
TcoU- We sec that - as we would expect - only the relative position is affected by the 
collision term. For this reason, it's useful to change coordinate to the centre of mass 
and the relative positions of the two particles. We write 

1 

R^2^n + f2) , f=fi-f2 

and similar for the momentum 

-P = Pi+P2 , P=^(pi-P2) 
And we can think of /2 = f2{R, f, P,P', t)- The distribution function will depend on the 

— * — * 

centre of mass variables R and P in some slow fashion, much as fi depends on position 



^ d_ 

dpi dp2 



f2^0 (2.32) 
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and momentum. In contrast, the dependence of /2 on the relative coordinates r and 
p is much faster - these vary over the short distance scale and can change on a time 
scale of order Tcow- 

Since the relative distributions in /2 vary much more quickly that /i, we'll assume 
that /2 reaches equilibrium and then feeds into the dynamics of /i. This means that, 
ignoring the slow variations in R and P, we will assume that df2ldt = and replace 
(2.32) with the equilibrium condition 

/ p_ dUjr) d\j,^^^ .2.33) 

\ m dr df dp J 

This is now in a form that allows us to start manipulating the collision integral on the 
right-hand-side of (2.30). We have 

dt y^^u J dfi dpi 



df 



, d-'r2d-% (pi - P2) • ^ (2.34) 

m 7|n-r-2|<d or 

where in the second line the extra term d/dp^ vanishes if we integrate by parts and, 
in the third line, we've used our equilibrium condition (2.33), with the limits on the 
integral in place to remind us that only the region r < d contributes to the collision 
integral. 

A Review of Scattering Cross Sections 

To complete the story, we still need to turn (2.34) into the coUision integral (2.24). 
But most of the work simply involves clarifying how the scattering rate <^(p,P2|Pi,P2) 
is defined for a given inter-atomic potential U{fi—f2). And, for this, we need to review 
the concept of the differential cross section. 

Let's think about the coUision between two particles. They start with momenta 
Pi — mvi and end with momenta = mt?/ with i — 1,2. Now let's pick a favourite, 
say particle 1. We'll sit in its rest frame and consider an onslaught of bombarding 
particles, each with velocity V2 — V1. This beam of incoming particles do not all hit our 
favourite boy at the same point. Instead, they come in randomly distributed over the 
plane perpendicular to V2 — Vi. The flux, /, of these incoming particles is the number 
hitting this plane per area per second, 

/ = y\V2 - Vi\ 
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Figure 4: The differential cross section. 



Now spend some time staring at Figure 4. There are a number of quantities defined 
in tliis picture. First, the impact parameter, b, is the distance from the asymptotic 
trajectory to the dotted, centre hne. We will use b and (p as polar coordinates to 
parameterize the plane perpendicular to the incoming particle. Next, the scattering 
angle, 6, is the angle by which the incoming particle is deflected. Finally, there are two 
solid angles, da and dQ, depicted in the figure. Geometrically, we see that they are 
given by 

da — bdbd(f) and = sin OdOdcf) 
The number of particles scattered into dfl in unit time is Ida. We usually write this as 

(2.35) 

where the differential cross section is defined as 

d{b^) 



I'^dVt^ Ibdbd(t) 
dU 



da 



dQ 



sin 9 



db 



de 



dcos9 



(2.36) 



You should think of this in the following way: for a fixed {v2 — Vi), there is a unique 
relationship between the impact parameter b and the scattering angle 6 and, for a given 
potential U{r), you need to figure this out to get \da/dQ\ as a function of 9. 
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Figure 5: Hard sphere scattering. 



Now we can compare this to the notation that we used earher in (2.17). There we 
talked about the rate of scattering into a small area (Pp'^d^p^ in momentum space. But 
this is the same thing as the differential cross-section. 



V2 



da 



dn 



dfl 



(2.37) 



(Note, if you're worried about the fact that d^p[d^p2 is a six-dimensional area while 
dD, is a two dimensional area, recall that conservation of energy and momenta provide 
four restrictions on the abihty of particles to scatter. These are imphcit on the left, 
but explicit on the right). 

An Example: Hard Spheres 

In Section 1.2, we modelled atoms as hard spheres of diameter d. Suppose that they 
scatter at impact parameter b < d. From the figure on the right, we see that the 
scattering angle is 6* = tt — 2q; where b = {d/2){l + sin a). So we have 

6^ = ^(l + sin^(|-i^)) = ^(2 + cos^) 
So, from (2.36), we find the differential cross-section 



da 



dfl 



The total cross-section is defined as 



ax 



16 

1 

27r / de sine^ 
Iq dil 



TT 



which provides a nice justification for the name because this is indeed the cross-sectional 
area of a sphere of radius d/2. 
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Figure 6: Two particle scattering. 



Almost Done 

With this refresher course on classical scattering, we can return to the collision integral 
(2.34) in the Boltzmann equation. 

V ^'^ / coll J\ri-f2\<d or 

We'll work in cylindrical polar coordinates shown in Figure 6. The direction parallel 
to V2 — vi is parameterized by x; the plane perpendicular is parameterised by (p and 
the impact parameter b. We've also shown the collision zone in this figure. Using the 
definitions (2.35) and (2.37), we have 

= j d^P2d%d% Uj{p[,p^\p,p2) [f2{x2) - f2{Xl)] 

It remains only to decide what form the two-particle distribution function /2 takes just 
before the collision at x = Xi and just after the collision at At this point we 

invoke the assumption of molecular chaos. Just before we enter the collision, we assume 
that the two particles are uncorrelated. Moreover, we assume that the two particles 
are once again uncorrelated by the time they leave the collision, albeit now with their 
new momenta 

f2{xi) = fi{f,pi;t)fi{f,p2;t) and /2(x2) = fi{r,p[;t)fi{f,p^;t) 
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Notice that all functions /i are evaluated at the same point r in space since we've 
assumed that the single particle distribution function is suitably coarse grained that it 
doesn't vary on scales of order d. With this final assumption, we get what we wanted: 
the collision integral is given by 



in agreement with (2.24). 
2.3 The H-Theorem 

The topics of thermodynamics and statistical mechanics are all to do with the equi- 
librium properties of systems. One of the key intuitive ideas that underpins their 

importance is that if you wait long enough, any system will eventually settle down to 
equilibrium. But how do wc know this? Moreover, it seems that it would be rather 
tricky to prove: settling down to equilibrium clearly involves an arrow of time that dis- 
tinguishes the future from the past. Yet the underlying classical mechanics is invariant 
under time reversal. 

The purpose of this section is to demonstrate that, within the framework of the 
Boltzmann equation, systems do indeed settle down to equilibrium. As we described 
above, wc have introduced an arrow of time into the Boltzmann equation. We didn't 
do this in any crude way like adding friction to the system. Instead, we merely assumed 
that particle velocities were uncorrelated before coUisions. That would seem to be a 
rather minor input but, as we will now show, it's enough to demonstrate the approach 
to equilibrium. 

Specifically, we will prove the "H-theorem", named after a quantity H introduced 
by Boltzmann. (H is not to be confused with the Hamiltonian. Boltzmann originally 
called this quantity something like a German £, but the letter was somehow lost in 
translation and the name H stuck). This quantity is 



This kind of expression is familiar from our first statistical mechanics course where we 
saw that the entropy S for a probability distribution p is S — — /c^plogp. In other 
words, this quantity H is simply 



The H-theorem, first proven by Boltzmann in 1872, is the statement that H always 
decreases with time. The entropy always increases. We will now prove this. 




d^P2d%d^P2 uj{p[,P2\p,P2) fi{r,p[)fi{r,P2) - fi{r,p)fi{'r,P2) 




S = -ksH 
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As in the derivation (2.4), when you're looking at the variation of expectation values 
you only care about the explicit time dependence, meaning 



dH 
It 



d^rd^p (log/i + 1) 



dh 
dt 



d^rd^p log/i 



dt 



where we can drop the +1 because J fi — N is unchanging, ensuring that J dfi/dt — 0. 
Using the Boltzmann equation (2.23), we have 

— = fdhd'p lo fJ^.^-l.^+ (^] ] 
dt J \df dp m df \ dt J 

But the first two terms in this expression both vanish. You can see this by integrating 
by parts twice, first moving the derivative away from /i and onto log/i, and then 
moving it back. We learn that the change in H is governed entirely by the collision 
terms 



dH 
'dt 



d^rd^p log/: 



dt 



= / d^rd^pid-^p2d-^p[d-^P2 uj{p[,p2\pi,p2) \ogfi{pi) 



coll 



X 



/i(p1)/i(P2)-/i(Pi)/i(P2) 



(2.38) 



where I've suppressed r and t arguments of /i to keep things looking vaguely reasonable 
I've also relabelled the integration variable p — > pi. At this stage, all momenta are 
integrated over so they are really nothing but dummy variables. Let's relabel 1 <-> 2 on 
the momenta. All the terms remain unchanged except the log. So we can also write 



dH 
It 



dhd^p,d'p2d%d% u;{p[,p^\p,,p2) log/i(p2) 



X 



/i(p1)/i(P2)-/i(Pi)/i(P2) 



(2.39) 



Adding (2.38) and (2.39), we have the more symmetric looking expression 
dH 1 



dt 



dhd'p,d'p2d'p[d% ^(p1,P2|Pi,P2) log [/i(pi) /l(p2)] 



X 



/i(p1)/i(P2)-/i(Pi)/i(P2) 



(2.40) 



Since all momenta are integrated over, we're allowed to just flip the dummy indices 
again. This time we swap p <-> p' in the above expression. But, using the symmetry 
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property (2.20), the scattering function remains unchanged^. We get 

^ = y d\d'pi(fp2d'p\£p'^ ^{p[J'2\Plj2) l0g[/l(pD/l(P2)] 

X 'fM)fM) - /i(Pi)/i(P2)] (2.41) 
Finally, we add (2.40) and (2.41) to get 

^ = y (frd'pid'p2(fp[(fp'2 UJ{P[,P2\P1:P2) X 

"log [fM) fM)] - log [/i(pi) Mp2)] ] [fM)fM) - /i(Pi)/i(P2)] (2.42) 

The bottom line of this expression is a function (log a; — logy)(x — y). It is positive for 
all values of x and y. Since the scattering rate is also positive, we have the proof of the 
H-theorem. 

dH ^ dS ^ 

— < <^ -r > 
dt ^ dt - 

And there we see the arrow of time seemingly emerging from time-invariant Hamiltonian 
mechanics! Clearly, this should be impossible, a point first made by Loschmidt soon 
after Boltzmann's original derivation. But, as we saw earlier, everything hinges on the 
assumption of molecular chaos (2.22). This was where we broke time-reversal symmetry, 
ultimately ensuring that entropy increases only in the future. Had we instead decided 
in (2.21) that the rate of scattering was proportional to /2 after the collision, again 
assuming /2 ~ /i/i then we would find that entropy always decreases as we move into 
the future. 

There is much discussion in the literature about the importance of the H-theorem and 
its relationship to the second law of thermodynamics. Notably, it is not particularly 
hard to construct states which violate the H-theorem by virtue of their failure to obey 
the assumption of molecular chaos. Nonetheless, these states still obey a suitable second 
law of thermodynamics^. 



^An aside: it's not actually necessary to assume (2.20) to make this step. We can get away with 
the weaker result 

j (fp[d^P2 L0{pi,P2\Pl,P2) = j (fp'id^P2 w(pi,P2|Pl,P2) 

which follows from unitarity of the scattering matrix. 

''This was first pointed out by E. T. Jaynes in the paper " Violation of Boltzmann's H Theorem in 
Real Gasef , published in Physical Review A, volume 4, number 2 (1971). 
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The H-theorem is not a strict inequality. For some distributions, the entropy remains 
unchanged. From (2.42), we see that these obey 



But this is simply the requirement of detailed balance (2.25). And, as we have seen al- 
ready, this is obeyed by any distribution satisfying the requirement of local equilibrium 
(2.29). 

2.4 A First Look at Hydrodynamics 

Hydrodynamics is what you get if you take thermodynamics and splash it. You know 
from your first course on statistical mechanics that, at the most coarse grained level, 
the equilibrium properties of any system are governed by the thermodynamics. In the 
same manner, low energy, long wavelength, excitations of any system are described by 
hydrodynamics. 

More precisely, hydrodynamics describes the dynamics of systems that are in local 
equilibrium, with parameters that vary slowly in space in time. As we will see, this 
means that the relevant dynamical variables are, in the simplest cases, 

• Density p{r,t) = mn{f,t) 

• Temperature r(r, t) 

• Velocity u{f, t) 

Our goal in this section is to understand why these are the relevant variables to describe 
the system and to derive the equations that govern their dynamics. 

2.4.1 Conserved Quantities 

We'll start by answering the first question: why are these the variables of interest? The 
answer is that these are quantities which don't relax back down to their equilibrium 
value in an atomic blink of an eye, but instead change on a much slower, domestic time 
scale. At heart, the reason for they have this property is that they are all associated 
to conserved quantities. Let's see why. 

Consider a general function A{f,p) over the single particle phase space. Because we 
live in real space instead of momentum space, the question of how things vary with 
r is more immediately interesting. For this reason, we integrate over momentum and 
define the average of a quantity v4(f, to be 



fM)fM)-fi{pi)fi{p2) 
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However, we've already got a name for the denominator in this expression: it is the 
number density of particles 



n{f,t) = J d'ph{f,p;t) (2.43) 

(As a check of the consistency of our notation, if you plug the local equilibrium dis- 
tribution (2.29) into this expression, then the n{r,t) on the left-hand-side equals the 
n{r,t) defined in (2.29)). So the average is 

(A(f, t)) = j d^p A{f,p)h{f,p; t) (2.44) 

It's worth making a couple of simple remarks. Firstly, this is different from the average 
that we defined earlier in (2.3) when discussion Liouville evolution. Here we're inte- 
grating only over momenta and the resulting average is a function of space. A related 
point is that we're at liberty to take functions which depend only on r (and not on p) 
in and out of the (•) brackets. So, for example, (nA) — n{A). 

We're interested in the how the average of A changes with time. We looked at this 
kind of question for Liouville evolution earlier in this section and found the answer 
(2.5). Now we want to ask the same question for the Boltzmann equation. Before we 
actually write down the answer, you can guess what it will look like: there will be a 
streaming term and a term due to the collision integral. Moreover, we know from our 
previous discussion that the term involving the collision integral will vary much faster 
than the streaming term. 

Since we're ultimately interested in quantities which vary slowly, this motivates look- 
ing at functions A which vanish when integrated against the collision integral. We will 
see shortly that the relevant criterion is 



We'd like to find quantities A which have this property for any distribution fi. Using 
our expression for the collision integral (2.23), we want 



d^Pid^P2d^p[d^p'2 uj{p[,p2\p,P2) A{f,pi) fi{f,p[)fi{f,P2) - fi{r,p)fi{f,P2 







This now looks rather similar to equation (2.38), just with the log / replaced by A. In- 
deed, we can follow the steps between (2.38) and (2.41), using the symmetry properties 
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of a;, to massage this into the form 



d^Pld^P2d%d% Ujip[,P2\Pl,P2) [/i(pO/i(P2) -/l(Pl)/l(P2) 

'A{f,p^) + A{f,p2) - A{f,pi) - A{f,pi) 



X 







Now it's clear that if we want this to vanish for all distributions, then A itself must 
have the property that it remains unchanged before and after the collision, 

A{f, p,) + A(f, P2) = A(f, p[) + A(f, p1) (2.45) 

Quantities which obey this are sometimes called collisional invariants. Of course, in 
the simplest situation we already know what they are: momentum (2.18) and energy 
(2.19) and, not forgetting, the trivial solution A — 1. We'll turn to each of these in 
turn shortly. But first let's derive an expression for the time evolution of any quantity 
obeying (2.45). 

Take the Boltzmann equation (2.23), multiply by a collisional invariant A{f,p) and 
integrate over J d^p. Because the collision term vanishes, we have 

/d»p.4(f,p)(| + |.A + /.|)/,(f.p,t) = 

where the external force is F = — W. We'll integrate the last term by parts (remem- 
bering that the force F can depend on position but not on momentum). We can't 
integrate the middle term by parts since we're not integrating over space, but nonethe- 
less, we'll also rewrite it. Finally, since A has no explicit time dependence, we can take 
it inside the time derivative. We have 



dt J df J m J m df 

Although this doesn't really look like an improvement, the advantage of writing it in 
this way is apparent when we remember our expression for the average (2.44). Using 
this notation, we can write the evolution of A as 

f) r) r)A r)A 

-{nA) + - . {nvA) - n{v ■ -) - n{F ■ -) ^ (2.46) 

where v = p/m. This is our master equation that tells us how any collisional invariant 
changes. The next step is to look at specific quantities. There are three and we'll take 
each in turn 
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Density 

Our first collisional invariant is the trivial one: A — 1. If we plug this into (2.46) we 
get the equation for the particle density n{f, t), 

9n d , ^. . 

(2.47) 

where the average velocity u of the particles is defined by 

u{t, t) = (v) 

Notice that, once again, our notation is consistent with earlier definitions: if we pick 
the local equilibrium distribution (2.29), the u{f,t) in (2.29) agrees with that defined 
above. The result (2.47) is the continuity equation, expressing the conservation of 
particle number. Notice, however, that this is not a closed expression for the particle 
density n: we need to know the velocity u as well. 

It's useful to give a couple of extra, trivial, definitions at this stage. First, although 
we won't use this notation, the continuity equation is sometimes written in terms of 

— * 

the current, J{r,t) = n{f,t) u{f,t). In what follows, we will often replace the particle 
density with the mass density. 



p(f, t) = mn{f, t) 



Momentum 



Out next collisional invariant is the momentum. We substitute A — mv into (2.46) to 
find 

d d 

— {mnui) + —{mnVjVi) - (nFi) = (2.48) 

We can play around with the middle term a little. We write 

{vjVi) = {{vj - Uj){vi - Ui)) + Ui{vj) + Uj{vi) - iiUj 

= {{Vj - Uj){Vi - Ui)) + UiUj 

We define a new object known as the pressure tensor, 

Pij = Pji = P((^j - Uj){Vi - Ui)) 

This tensor is computing the flux of i-momentum in the j-direction. It's worth pausing 
to see why this is related to pressure. Clearly, the exact form of Pij depends on the 
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distribution of particles. But, we can evaluate the pressure tensor on the equilibrium, 
Maxwell-Boltzmann distribution (2.28). The calculation boils down to the same one 
you did in your first Statistical Physics course to compute equipartition: you find 



Pij = nkBT5ij 



(2.49) 



which, by the ideal gas law, is proportional to the pressure of the gas. Using this 
definition - together with the continuity equation (2.47) - we can write (2.48) as 



This is the equation which captures momentum conservation in our system. Indeed, it 
has a simple interpretation in terms of Newton's second law. The left-hand-side is the 
acceleration of an clement of fluid. The combination of derivatives is sometimes called 
the material derivative, 



It captures the rate of change of a quantity as seen by an observer swept along the 
streamline of the fluid. The right-hand side of (2.50) includes both the external force 
F and an additional term involving the internal pressure of the fluid. As we will see 
later, ultimately viscous terms will also come from here. 

Note that, once again, the equation (2.50) docs not provide a closed equation for 
the velocity u. You now need to know the pressure tensor P^j which depends on the 
particular distribution. 

Kinetic Energy 

Our final coUisional invariant is the kinetic energy of the particles. However, rather 
than take the absolute kinetic energy, it is slightly easier if we work with the relative 
kinetic energy. 



2 ' ' 

If we substitute this into the master equation^ (2.46), the term involving the force 
vanishes (because (vj — Ui) — 0). However, the term that involves dE/dri is not zero 

^There is actually a subtlety here. In deriving the master equation (2.46), we assumed that A has 

no explicit time dependence, but the A defined above does have exphcit time dependence through 
u{f,t). Nonetheless, you can check that (2.46) still holds, essentially because the extra term that you 
get is ~ {{v — u) ■ du/dt) — {v — u) ■ du/dt = 0. 




(2.50) 
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because the average velocity u depends on r. We have 

liiPi^-^'^f) + li:iP^^i^-^^f)-Pi-^^^i^-^'^f) = (2.52) 

At this point, we define the temperature, T{f, t) of our non-equihbrium system. To do 
so, we fall back on the idea of equipartition and write 

^mkBT{f,t) = ^{{v-u{f,t)f) (2.53) 

This coincides with our familiar definition of temperature for a system in local equilib- 
rium (2.29), but now extends this to a system that is out of equilibrium. Note that the 
temperature is a close relative of the pressure tensor, TrP = SmpkBT. 

We also define a new quantity, the heat flux, 

Qi = ^mpdvi - Ui) {v - uf) (2.54) 

(This actually difi'ers by an overall factor of m from the definition of q that we made in 
Section 1. This has the advantage of making the formulae we're about to derive a little 
cleaner). The utility of both of these definitions becomes apparent if we play around 
with the middle term in (2.52). We can write 

]^mp{vi{v - uf) = ^mp{{vi - Ui) {v - uf) + ]^mpUi{{v - uf) 
= qi + ^pUiksT 

Invoking the definition of the pressure tensor (2.49), we can now rewrite (2.52) as 
3 9 , , d f 3 , _\ „ duj 

Because Pij = Pji, we can replace duj/dri in the last term with the symmetric tensor 
known as the rate of strain (and I promise this is the last new definition for a while!) 

Finally, with a little help from the continuity equation (2.47), our expression for the 
conservation of energy becomes 
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It's been a bit of a slog, but finally we have three equations describing how the particle 
density n (2.47), the velocity u (2.50) and the temperature T (2.56) change with time. 
It's worth stressing that these equations hold for any distribution /i. However, the 
set of equations are not closed. The equation for n depends on u\ the equation for u 
depends on P^j and the equation for T (which is related to the trace of P^) depends 
on a new quantity q. And to determine any of these, we need to solve the Boltzmann 
equation and compute the distribution /i. But the Boltzmann equation is hard! How 
to do this? 

2.4.2 Ideal Fluids 

We start by simply guessing a form of the distribution function /i(r, p; t). We know that 
the collision term in the Boltzmann equation induces a fast relaxation to equilibrium, 
so if we're looking for a slowly varying solution a good guess is to take a distribution 
for which {dfi/dt)co\\ — 0. But we've already met distribution functions that obey this 
condition in (2.29): they are those describing local equilibrium. Therefore, our first 
guess for the distribution, which we write as ff'\ is local equilibrium 

/;«'(r-p;0 =n(f.O (^_J^)"'exp (-^^ t)]^) (2.57) 

where p = mv. In general, this distribution is not a solution to the Boltzmann equation 
since it does not vanish on the streaming terms. Nonetheless, we will take it as our 
first approximation to the true solution and later see what we're missing. 

The distribution is normalized so that the number density and temperature defined 
in (2.43) and (2.53) respectively coincide with n{f,t) and T{f,t) in (2.29). But we can 
also use the distribution to compute and q. We have 

Pij = kBn{f, t)T{r, t) Sij = P{f, t) 5ij (2.58) 

and q = 0. We can substitute these expressions into our three conservation laws. The 
continuity equation (2.47) remains unchanged. Written in terms for p = mn, it reads 

d d \ dui 

Meanwhile, the equation (2.50) governing the velocity fiow becomes the Euler equation 
describing fiuid motion 

d_ d\ ldP_ _F\ 
dt ^dfj) ^ pdvi m 
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and the final equation (2.56) describing the flow of heat reduces to 




(2.61) 



These set of equations describe the motion of an ideal fluid. While they are a good 
starting point for describing many properties of fluid mechanics, there is one thing that 
they are missing: dissipation. There is no irreversibility sown into these equations, no 
mechanism for the fluid to return to equilibrium. 

We may have anticipated that these equations lack dissipation. Their starting point 
was the local equilibrium distribution (2.57) and we saw earlier that for such distribu- 
tions Boltzmann's H-function does not decrease; the entropy does not increase. In fact, 
we can also show this statement directly from the equations above. We can combine 



which tells us that the quantity pT~^/^ is constant along streamlines. But this is the 
requirement that motion along streamlines is adiabatic, not increasing the entropy. To 
see that this is the case, you need to go back to your earlier statistical mechanics or 
thermodynamics course^. The usual statement is that for an ideal gas, an adiabatic 
transformation leaves VT'^^^ constant. Here we're working with the density p — mN/V 
and this becomes pT~'^/'^ is constant. Note, however, that in the present context p and 
T are not numbers, but functions of space and time: we are now talking about a local 
adiabatic change. 

Sound Waves 

It is also simple to show explicitly that one can set up motion in the ideal fluid that 
doesn't relax back down to equilibrium. We start with a fluid at rest, setting u — Q 
and p — p and T — T, with both p and T constant. We now splash it (gently). That 
means that we perturb the system and linearise the resulting equations. We'll analyse 
these perturbations in Fourier modes and write 



Furthermore, we'll look for a particular kind of perturbation in which the fluid motion 
is parallel to the perturbation. In other words, we're looking for a longitudinal wave 



^See, for example, the discussion of the Carnot cycle in Section 4 of the lecture notes on Statistical 
Physics: http://www.damtp.cani.ac.uk/user/tong/statphys.htnil 



(2.59) and (2.60) to flnd 




p(f , t)^p + 5p e-'^'^'-''-^ and T(f, t)^f + 5T e''^'^'-''-^ 



(2.62) 




(2.63) 
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The hnearised versions of (2.59), (2.60) and (2.61) then read 



So- 

\k\ 


p5u 


^ A., 
-^OU = 

\k\ 


^bT kBj- 

op H 01 

mp m 


^ST = 

\k\ 


IfSu 



There is one solution to these equations with zero frequency, u = 0. These have Su = 
while 6p = —p and 6T = T. (Note that this notation hides a small e. It really means 
that Sp = —ep and 5T = eT. Because the equations are linear and homogeneous, you 
can take any e you like but, since we're looking at small perturbations, it should be 
small). This solution has the property that P — mnkBT is constant. But since, in 
the absence of an external force, pressure is the only driving term in (2.60), the fluid 
remains at rest, which is why 5u — Qior this solution. 

Two further solutions to these equations both have 5p = p, 6T = |T and 6u = u!/\k\ 
with the dispersion relation 

/ 5 T 

(J — ±Vs\k\ with Vs — \ — — (2.64) 

V 3m 

These are sound waves, the propagating version of the adiabatic change that we saw 
above: the combination pT~^/^ is left unchanged by the compression and expansion of 
the fluid. The quantity Vs is the speed of sound. 

2.5 Transport with Collisions 

While it's nice to have derived some simple equations describing fluid mechanics, as 
we've seen they're missing dissipation. And, since the purported goal of these lectures 
is to understand how systems relax back to equilibrium, we should try to see what 
we've missed. 

In fact, it's clear what we've missed. Our flrst guess for the distribution function was 
local equilibrium 

fr(f,m ^ n(r,t) (^--J^)^^xp (-^^ S(f,t)f) (2.65) 

We chose this on the grounds that it gives a vanishing contribution to the collision 
integral. But we never checked whether it actually solves the streaming terms in the 
Boltzmann equation. And, as we will now show, it doesn't. 
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Using the definition of the Poisson bracket and the one-particle Hamiltonian Hi 
(2.11), we have 



(0) 



{Huff} 



(0) 



+ F 



(0) 



(0) 



dt dt dp df 

Now the dependence on p = mv in local equilibrium is easy: it is simply 



(0) 



(0) 



Meanwhile all r dependence and t dependence of f[^^ lies in the functions n(r, t), T(f, t) 
and u{f,t). From (2.65) we have 
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Using all these relations, we have 
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knT 



f?' (2.66) 



where we've introduced the notation Df which differs from the material derivative 
in that it depends on the velocity v rather than the average velocity u, 

d 



d_ ^ d_ 
dt df 



Dt + {v-u) 



df 



Now our first attempt at deriving hydrodynamics gave us three equations describing 
how n (2.59), u (2.60) and T (2.61) change with time. We substitute these into (2.66). 
You'll need a couple of lines of algebra, cancelling some terms, using the relationship 
P — nksT and the definition of Uij in (2.55), but it's not hard to show that we 
ultimately get 

Off' 



dt 



^ ( -4^{v-uf-^] {v-u)-VT 
T \ 2kBT^ ' 2 ' ^ ' 



(2.67) 
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And there's no reason that the right-hand-side is zero. So, unsurprisingly, /} does 
not solve the Boltzmann equation. However, the remaining term depends on VT and 
du/dr which means that we if we stick to long wavelength variations in the temperature 
and velocity then we almost have a solution. We need only add a httle extra something 
to the distribution 



Let's see how this changes things. 

2.5.1 Relaxation Time Approximation 

The correction term, will contribute to the collision integral (2.24). Dropping the 
r argument for clarity, we have 



where, in the second line, we have used the fact that /{ vanishes in the collision 
integral and ignored quadratic terms ~ . The resulting collision integral is a hnear 
function of 5fi. But it's still kind of a mess and not easy to play with. 

At this point, there is a proper way to proceed. This involves first taking more care 
in the expansion of 5fi (using what is known as the Chapman- Enskog expansion) and 
then treating the linear operator above correctly. However, there is a much easier way 
to make progress: we just replace the collision integral with another, much simpler 
function, that captures much of the relevant physics. We take 



where r is the relaxation time which, as we've already seen, governs the rate of change 
of /i. In general, r could be momentum dependent. Here we'll simply take it to be a 
constant. 

The choice of operator (2.69) is called the relaxation time approximation. (Sometimes 
it is referred to as the Bhatnagar-Gross-Krook operator). It's most certainly not exact. 




(2.68) 




'p2rfVlrfV2 ^(p1,p1|Pi,P2) [/i(p1)/i(p1) - /l(Pl)/l(P2)] 



j d'p^dy,d'p',u{p[,p^\p,J,) [A'\p[)Sh{p^ + 5f{p[)f?\p^ 

-ff\n)5HP2) - 6f{pi)fi'\p2) 




(2.69) 
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In fact, it's a rather cheap approximation. But it will give us a good intuition for what's 
going on. With this replacement, the Boltzmann equation becomes 



{H,J, — 

But, since Sfi <^ ff'\ we can ignore 6fi on the left-hand-side. Then, using (2.67), we 
have a simple expression for the extra contribution to the distribution function 



Sfi = -r 



m 



T V 2kBT 

m 



[v-uf- 



[V — u) 



dT 

dr 

1.^ 



/r (2.70) 



We can now use this small correction to the distribution to revisit some of the transport 
properties that we saw in Section 1. 



2.5.2 Thermal Conductivity Revisited 

Let's start by computing the heat flux 

1 



(li 



mp{{vi - Ui) (v-u) ) 



(2.71) 



using the corrected distribution (2.68). We've already seen that the local equihbrium 
distribution /}°^ gave q — 0, so the only contribution comes from 5fi. Moreover, only 
the first term in (2.70) contributes to (2.71). (The other is an odd function and vanishes 
when we do the integral) . We have 

q= -kVT 

This is the same phenomenological law that we met in (1.12). The coefficient k is the 
thermal conductivity and is given by 



K 



mrp 
mrp 



m 2 5' 

2fc^("-") -2 



(0) 



In the second line, we've replaced all {v — u) factors with v by performing a {f- 
dependent) shift of the integration variable. The subscript {■)q means that these aver- 
ages are to be taken in the local Maxwell-Boltzmann distribution f^'^ with u — Q. These 
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integrals are simple to perform. We have {v'^)o = IS/c^T^/m^ and {v^)o = 105k%T^/m^, 
giving 

K — -rnk'iiT 

The factor oi 5/2 here has followed us throughout the calculation. The reason for its 
presence is that its the specific heat at constant pressure, Cp — ^ks- 

This result is parameterically the same that we found earlier in (1.13). (Although 
you have to be a httle careful to check this because, as we mentioned after (2.54), 
the definition of heat flux differs and, correspondingly, k, differs by a factor of m. 
Moreover, the current formula is written in terms of slightly different variables. To 
make the comparison, you should rewrite the scattering time as r ~ l/man^JJv^ , 
where a is the total cross-section and (f^) ~ T/m by equipartition). The coefficient 
differs from our earlier derivation, but it's not really to be trusted here, not least 
because the only definition of r that we have is in the implementation of the relaxation 
time approximation. 

We can also see how the equation (2.56) governing the flow of temperature is related 
to the more simplistic heat flow equation that we introduced in (1.14). For this we 
need to assume both a static fluid -u = and also that we can neglect changes in the 
thermal conductivity, dn/df'^ 0. Then equation (2.56) reduces to the heat equation 

, dT 2 

2.5.3 Viscosity Revisited 

Let's now look at the shear viscosity. From our discussion in Section 1, we know that 
the relevant experimental set-up is a fluid with a velocity gradient, du^/dz ^ 0. The 
shear viscosity is associated to the flux of x-momentum in the ^-direction. But this is 
precisely what is computed by the off-diagonal component of the pressure tensor, 

Pxz = p{{Vx - U:c){Vz - Uz)) 

We've already seen that the local equilibrium distribution gives a diagonal pressure 
tensor (2.58), corresponding to vanishing viscosity. What happens if we use the cor- 
rected distribution (2.68). Now only the second term in (2.70) contributes (since the 
flrst term is gives an odd function of {v — u)). We write 



(2.72) 



where the extra term Iljj is called the stress tensor and is given by 



LLij — , ^ Ukl 
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kl 



d p {vj - Uj){vi - Ui) {vk - ui){vk - ui) - - {v-u) 6ki /{ 



{ViVjVkVi)o - -Skl{ViVjV^)o 



p(0) 



Before we compute 11^, note that it is a traceless tensor. This is because the first 
term above becomes {v^VkVi)o — Sjk{v^VxVx)o which is easily calculated to be {v^v'^)q — 
bk\T'^ /m? = \(v'^)q. Moreover, Hij depends linearly on the tensor Uij. These two facts 
mean that 11^ must be of the form 



(2.73) 



In particular, if we set up a fluid gradient with dux/dz ^ 0, we have 

dux 



-Tj- 



dz 



which tells us that we should identify 77 with the shear viscosity. To compute it, we 
return to a general velocity profile which, from (2.73), gives 



mrp 
mrp 
2mrp 



kl 



{VxVzVkVl)o - -5kl{VxVzV'^ 



{Uxz + U^x){VxVzVxVz)o 
Uxz{v'^)o 



IbksT 

Comparing to (2.73), we get an expression for the coefficient r], 

r) — nksTr 

Once again, this differs from our earlier more naive analysis (1.11) only in the overall 
numerical coefficient. And, once again, this coefficient is not really trustworthy due to 
our reliance on the relaxation time approximation. 

The scattering time r occurs in both the thermal conductivity and the viscosity. Tak- 
ing the ratio of the two, we can construct a dimensionless number which characterises 
our system. This is called the Prandtl number, 

Pr = M 

K 
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With Cp the specific heat at constant pressure which takes the value Cp = 5A;b/2 for a 
monatomic gas. Our calculations above give a Prandtl number Pr = 1. Experimental 
data for monatomic gases shows a range of Prandtl numbers, hovering around Pr 2/3. 
The reason for the discrepancy lies in the use of the relaxation time approximation. A 
more direct treatment of the collision integral, thought of as a linear operator acting 
on Sfi, gives the result Pr = 2/3, in much better agreement with the data^. 

2.6 A Second Look: The Navier-Stokes Equation 

To end our discussion of kinetic theory, we put together our set of equations governing 
the conservation of density, momentum and energy with the corrected distribution 
function. The equation of motion for density fluctuations doe not change: it remains, 

1^ + V • (pu) = (2.74) 

Meanwhile the equation for momentum (2.50) now has an extra contribution from the 
stress tensor contribution (2.72). Moreover, we typically assume that, to leading order, 
variations in the viscosity can be neglected: Vr] ~ 0. Written in vector notation rather 
than index notation, the resulting equation is 

^ + u-v]u^-- -VP + ^V'lT + :^V(V • u) (2.75) 
at J m p p 3p 

This is the Navier-Stokes equation. Finally, we have the heat conduction equation. We 
again drop some terms on the grounds that they are small. This time, we set V«; ~ 
and UijUij 0; both are small at the order we are working to. We're left with 

/d \ 2 , 2m 

We can again look at fluctuations of these equations about a static fluid with p = p, 
T — T and u — 0. Longitudinal fluctuations (2.62) and (2.63) now give rise to the 
linearised equations of motion, 

LuSp = p\k\6u 

.5u='^\k\Sp+^-^\k\5T-^^-^Su 
mp m 6p 

2 - ^ 2kIA-I^ 
uST^ -f\k\5u-i—^ST 

3 'iksp 



^You can read about this improved calculations in the lecture notes by Arovas. A link is given on 
the course website. 
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Notice that terms involving transport coefficients rj and k each come with a factor of 
i; this is a sign that they will give rise to dissipation. To compute the frequencies of 
the different modes, it's best to think of this as an eigenvalue problem for u!/\k\; the 
coefficients of the various terms on the right-hand-side define a 3 x 3 matrix M, with 

detM = — ' '_ and Tr M ^ -i { -t] + -— ] ^ 
3 mp \3 oKb J p 

The product of the three eigenvalues is equal to det M. We know that for the ideal 
fluid, the eigenvalues are zero and u — ±Vs\k\ where Vs is the sound speed computed 
in (2.64). Let's ffist look at the eigenvalue that was zero, corresponding to fluctuations 
of constant pressure. Working to leading order in k and rj, we must have 

-v^A k\^ uj = det M ^ uj = - — -^\k\^ 

5 ksp 

The purely imaginary frequency is telling us that these modes are damped. The uj ~ 
i\k\'^ is characteristic of diffusive behaviour. 

The remaining two modes are related to the sound waves. These too will gain a 
dispersive contribution, now with 

UJ — ■^Vs\k\ — i'-f (2.76) 
Using the fact that the sum of the eigenvalues is equal to the trace, we find 

The fluctuations above are all longitudinal. There are also two shear modes, whose 
fluctuations are in a direction perpendicular to the velocity. It is simple to check that 

— * 

the linearised equations are solved hy 5p — 5T — Q and 6u- k, with the frequency given 

by 

7]\k\^ 
UJ — —I 



P 

Once again, we see that these modes behave diffusively. 
Navier Stokes Equation and Liquids 

Our derivation of the Navier-Stokes equation relied on the dilute gas approximation. 
However, the equation is more general than that. Indeed, it can be thought of as the 
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most general expression in a derivative expansion for momentum transport (subject to 
various requirements). In fact, there is one extra parameter that we could include: 

p{^^ + u-V^u = (^-VP + riV^u + (| + c)v(V-m) 

where C, is the hulk viscosity which vanished in our derivation above. Although the 
equation above governs transport in liquids, we should stress that first-principles com- 
putations of the viscosity (and also thermal conductivity) that we saw previously only 
hold in the dilute gas approximation. 
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